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ABSTRACT 

We describe an implicit prediction of the accretion disk models constructed 
by Wardle and Konigl (1990) for the circumnuclear disk (CND) of gas and dust 
near the Galactic center: supersonic ambipolar diffusion, an essential dynamical 
ingredient of the Wardle-Konigl disks, will cause the alignment of dust grains due 
to a process described by Roberge, Hanany, & Messinger (1995). We calculate 
synthetic maps of the polarized thermal emission which would be caused by 
ambipolar alignment in the preferred Wardle-Konigl model. Our maps are in 
reasonable agreement with 100 /iin polarimetry of the CND if we assume that 
the grains have shapes similar to those of grains in nearby molecular clouds and 
that the CND contains a disordered magnetic field in energy equipartition with 
its ordered field. 

Subject headings: dust — Galaxy: center — infrared: ISM: continuum — ISM: 
magnetic fields — MHD — polarization 
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1. INTRODUCTION 

The Galactic circumnuclear disk (CND) is a weakly-ionized ring of gas and dust which 
appears in various types of emission at galactocentric radii R = 1-10 pc (see Morris & 
Serabyn 1996 and Davidson 1996 for recent reviews). The origin and structure of the CND 
magnetic field, B, are particularly intriguing: mid- and far-infrared polarimetry ([Aitken et 



|al- 1986| ; |Werner et al. 1988| ; Hildebrand et al. 1990, 1993 [H90, H93]) show that B lies 
primarily in the plane of the disk and is roughly orthogonal to the poloidal field observed 
on larger scales. Werner et al. (1988) pointed out that a differentially-rotating disk would 
deform an initially poloidal field into a predominantly toroidal configuration. Subsequently, 
Wardle and Konigl (1990, [WK]) described self-similar hydromagnetic accretion disk models 
of the CND; with the addition of a small nonaxisymmetric distortion, the magnetic field 
predicted by their preferred model, gc2, is in reasonable agreement with the field inferred 
from far- infrared polarimetry (H93). 

In the WK models, the magnetic torques which remove angular momentum from the 
gas are mediated by supersonic ambipolar diffusion. Coincidentally, Roberge, Hanany, & 
Messinger (1995, [RHM]) showed that ambipolar diffusion causes grain alignment: the par- 
tial coupling of charged dust grains to the drifting field lines produces gas-grain streaming, 
which causes alignment of the grains via Gold's mechanism (|Gold 19"5*2|; |Purcell 1969|; [Pur 



|cell fc Spitzer 1971| ; [Lazarian 1994| ; RHM). The large magnitude of the "intrinsic" 100 /im 



polarization from grains aligned by supersonic ambipolar diffusion (e.g., 20% for oblate sil- 
icate grains with 2:1 axis ratios viewed with the most favorable magnetic field geometry; 
see RHM) leads us to conjecture that the ambipolar mechanism is the predominant cause of 
grain alignment in the CND. In §§2-3, we test this conjecture by calculating the polarized 
thermal emission that would be caused by ambipolar alignment in model gc2 and comparing 
our results to the 100 fim polarimetry of H90 and H93. 
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2. CALCULATIONS 

The largest optical depth through the CND is < 0.1 at 100/im; consequently the Stokes 
parameters, J, Q, and U, are well described by the transfer equations for optically thin 
emission, 

^L = n d B u (T d ) (C^ + CV) (1) 

^ = n d B u (T d ) {C x , -C y ,)F cos 2^ (2) 
as 

^ = n d B„(T d ) (0,, - C y .) F sin 2^ (3) 

(WK90), where s is distance along the line of sight. We integrated equations (|T|)-(0) to obtain 
the polarization magnitude, P = VQ 2 + U 2 / 1, and position angleB, 6 = | arctan(?7/Q), 
using a model of the position-dependent dust density, n<j, and temperature, derived by 
Davidson et al. (1992) from 50 and 90 «m photometry of the Galactic center. Our results 
are only weakly sensitive to uncertainties in and because the CND is optically and 
geometrically thin and the FIR emission is in the Ray leigh- Jeans limit. The other variables 
in equations (|])-(§) are determined by the alignment properties of the grains as follows. 

The quantities C x > and C y i are respectively the dust absorption cross sections for light 
polarized along the x and y' directions, where y 1 is parallel to the projection of the local 
magnetic field onto the plane of the sky, z points toward the observer, x = y Xz , and ip is 
the angle between y and north, measured counterclockwise from north. To compute C x > and 
Cy>, we modeled the grains as oblate spheroids with semiaxes a parallel to the symmetry axis 
and b (> a) perpendicular to the symmetry axis. We assumed that the angular momentum, 
J, of a rotating spheroid is aligned perfectly with its symmetry axis by the Barnett effect.i 



x With 9 = if the electric vector is polarized along the east-west direction. 
2 That is, we neglected the disalignment of J by thermal fluctuations ( [Lazarian 1994j ). 
This is a good approximation for the conditions of interest here (Lazarian & Roberge 1996). 
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In this approximation, 

C x , = C avs + l -Qj (C ± - C||) (4) 

and 

C y/ = C avg + \Qj (l - 3 cos 2 C) (C ± - C||) (5) 



( [Lee fc Draine 1985 ), where C\\ and C± are respectively the cross sections for light polarized 



parallel and perpendicular to the symmetry axis, C avg = (2C± + C||) /3, and ( is the angle 
between B and the plane of the sky. For a given a/b ratio, we computed C\\ and C± in 
the electric dipole approximation using Draine's (1987) dielectric function for astronomical 
silicates. Notice that P and 9 depend only on ratios of the grain cross sections, which are 
independent of grain size in the electric dipole limit. Thus, apart from its weak influence 
on the grain dynamics (discussed below), the value of b does not affect our calculations. We 
assumed for simplicity that the grain shape is independent of location and let a/b be a free 
parameter. 

The angles ( and ip are functions of the position-dependent disk magnetic field. We 
computed £ and ip from a numerical tabulation of model gc2 kindly supplied by Mark Wardle. 
The model is self-similar and to determine the hydro variables uniquely it is also necessary 
to specify the magnetic field, -Bo, the density, nno, and the radial (ivo) and azimuthal (u^o) 
velocity components of the neutral gas at some fiducial radius, r . We set r = 1 pc and 
chose B = 2mG, consistent with H I and OH Zeeman splitting measurements of the CND 
(reviewed by Morris & Serabyn 1996). We set nno = 10 5 cm~ 3 , near the lower end of the 
range of estimates for the density ( |Jackson et al. 1993 ), and assumed that v^q = 110 km s _1 



and v r = —19 km s 1 , consistent with kinematic models of the CND (Gusten et al. 1987 
Jackson et al. 1993| ). 



The alignment of J with respect to B is characterized by the quantity 

3 



cos 2 /3)-i 



(6) 
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where (3 is the angle between J and B and angle brackets denote the average for all grains. 
We calculated Qj from the theory of ambipolar alignment (RHM).i For supersonic drift 
speeds, Qj depends only on a/b and Ud/^th, where t>d is the gas-grain drift speed and v t t = 
^/2kT g /m is the gas thermal speed (RHM). We computed v& at each point using model gc2 
to predict the local hydro variables and a prescription for the grain dynamics (Draine 1980) 
which includes Lorentz and gas drag forces on the grains. For the purposes of calculating the 
grain dynamics only, we represented the grains as spheres with equivalent radii r eq = (ab 2 ) 1 ^ 3 
and arbitrarily set b = 0.1 /im. Changing b by a factor of 2 would have virtually no effect 
on the calculated polarizations. The WK models do not predict the gas temperature; we set 
T g = 300 K, consistent with observations of various molecular lines ([Harris et al. 1985|) . 

The factor F in equations represents the possible effects of a random magnetic 

field component. We set the total field to B^g+SB, where -B avg is the ordered field predicted 
by model gc2 and 6B is a random component. We assumed that SB is perpendicular to B avg 
(as would be the case if 8B is due to Alfvenic turbulence) and that 8B has an axisymmetric, 
Gaussian distribution of amplitudes (see |Myers fc Goodman 1991] ). For this model, it is 



straightforward to show that 

F(£)= » r^^-I, (7) 
V27T4 jo 1 + x z 2 

where the "turbulence parameter," £, is the ratio of the energies in the ordered and disordered 
field components. We assumed for simplicity that this ratio is independent of position and 
let £ be a free parameter. Note that F decreases monotonically from F(0) = 1 (no disordered 
field) to F(l) ~ 0.5 (disordered field in equipartition with the ordered field). 



including Davis-Greenstein alignment would increase Qj, and hence the polarizations, 
by an uncertain factor that depends on poorly-known magnetic properties of the grains. The 
increase would be £ 40% for grains composed of ordinary paramagnetic substances but could 
be larger for superparamagnetic grains. 
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3. Results 

We constructed synthetic maps of P t hr, the 100 /im polarization predicted by ambipolar 
alignment, for 3,136 parameter combinations sampled uniformly on the intervals < a/b < 1 
and < £ < 2. For each (a/6, £) combination, we calculated the Stokes parameters on a grid 
of sightlines uniformly spaced by 1" in right ascension and declination. Each map of/, Q, and 
U was then convolved with a 45" (FWHM) Gaussian beam to model the observed resolution 
(H93) and maps of P t hr were computed from the convolved Stokes parameters. Optimal 
values of a/b and £ were obtained by performing a \ 2 fit to P bs> the 100 /im polarization 
observed along N = 23 lines of sight (H90, H93). Our fits omitted 7 observed sightlines for 
which the FIR emission is not associated with the CND. We did not attempt to fit 6, which 
is a good diagnostic of the magnetic field geometry but not of the alignment mechanism. 

Strictly speaking, the value of P Q bs we calculate for each sightline is the mean polariza- 
tion for a hypothetical ensemble of disks with identical distribution functions for SB. That 
is, discrepancies between P thr and P Q b s are due not only to uncertainties in our theoretical 
model but also to random fluctuations in B. To include the effects of random fluctuations 
on x 2 , we assumed that the polarization observed along the ith sightline has dispersion 
°f = °obs,i + °frb5 where cr b s ,i is the observational uncertainty quoted by H93 and <T tr b rep- 
resents the uncertainty due to fluctuations in B. We assumed for simplicity that cr tr b is 
the same for every sightline and performed the y 2 fits independently for different choices of 
Ctrb- We found that x 2 P er degree of freedom is < 1 in some subset of parameter space if 
o'trb > 0.6%. That is, we are able to obtain a reasonable fit to the data if we assume that 
the random fluctuations in P Q b s are about half a percent E 

4 This logic probably overestimates cr tr b, since some of the errors that are ascribed to 
fluctuations in P obs are undoubtedly due to idealizations in our model. For example, H93 
pointed out that the WK models predict polarizations that are symmetric in the northwest 
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Figure 1 is a contour plot of \ 2 P er degree of freedom (x 2 /^) f° r the case Ctrb = 0.6%. 
Evidently the observations are consistent with a continuum of models with different com- 
binations of a/6 and £; this is due to the fact that a small increase in a/6 (which tends to 
reduce Pthr) can always be compensated by a decrease in £ (which tends to increase Pthr)- 
Thus, the contours of constant x 2 jv slope downward with increasing a/6. It is interesting 
to note that the parameters allowed by our fits are consistent with independent constraints: 
For example, observations show nearby molecular clouds have £ ~ 1 (Jones 1989; Myers & 
Goodman 1991). Also, Zeeman splitting observations of the 21cm HI flSchwarz fc Lazcnby 



1990| ) and 18 cm OH ( [Killeen et al. 1993|) lines toward the Galactic center imply that the 
magnetic energy density of the CND is comparable to the turbulent energy density inferred 
from atomic and molecular linewidths flGtisten et al. 1987| ; | Jackson et al. 1993| ). These 



observations are consistent with the presence of nonlinear Alfven waves in the CND, for 
which £ « 1. Nonlinear Alfven waves are predicted by some models of magnetic braking 
(Mouschovias Sz Morton 1985). Similarly, photopolarimetry of the 9.7/xm silicate resonance 
implies that the aligned grains in nearby molecular clouds are oblate with a/6 rs 2/3 (Hilde- 
brand & Dragovan 1995). A model with a/6 = 2/3 and £ = 1 (indicated by the filled circle 
in Fig. 1) is compared with the observations in Figure 2. For this model, x 2 l v = 0.97 and 
the mean discrepancy between P t h r and P b s is —0.2%. 

Our results are insensitive to the precise values of the hydrodynamic variables; for 
example, the polarizations we calculate would be virtually identical if we arbitrarily increased 
or decreased hh or B by factors of ~ 3. This robustness is due to the fact that, for the highly 
supersonic drift speeds in the WK models (t>d ^ 10 for gc2), the efficiency of ambipolar 



and southeast quadrants. However, the observations clearly show departures from symmetry 
in the southeast quadrant; not surprisingly, these points contribute about twice as much on 
average to x 2 as the res t of the data. 
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alignment saturates at a value which depends only on the grain shape (see RHM, Fig. 12). 
However, our results are sensitive to the magnetic field geometry. For example, we attempted 
to fit the observations by arbitrarily replacing the magnetic field predicted by model gc2 with 
a purely toroidal field of the same magnitude at each point; we were unable to find solutions 
with x 2 / u < 2 for these models. It will be interesting to see whether one can model the 
CND observations successfully with more detailed hydro models (e.g., Wardle & Konigl 1993) 
which calculate the vertical structure of the disk. 

Previous models of the CND magnetic field geometry (WK, H90, H93) have accounted 
successfully for the observed polarization position angles by assuming that the efficiency 
of grain alignment is independent of position. The saturation of ambipolar alignment at 
large drift speeds provides a natural justification for this assumption. In contrast, it is 
difficult to reconcile the uniform efficiency required by the observations with other alignment 
mechanisms. If the grains are superparamagnetic, then the efficiency of Davis-Greenstein 
alignment would be saturated but the saturation efficiency would be very sensitive to the 
local dust-to-gas temperature ratio: independent calculations on DG alignment (Lazarian 
1995) show that factor of 2 variations in T d /T g would cause factor of ~ 3 variations in 
(cos 2 f3). The efficiency of Purcell's mechanism (Purcell 1979) is insensitive to T d /T g but 
sensitive to the timescale for changes in the grain surface properties; it is difficult to see why 
the latter should be uniform throughout the CND and have just the right value to reproduce 
the observations. We conclude that the degree and uniformity of the alignment required by 
the observations are natural consequences of ambipolar alignment but difficult to explain 
with other mechanisms. 
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4. SUMMARY 

1. We have calculated the 100 /xm linear polarization from the CND using a model that 
attributes the grain alignment to ambipolar diffusion in the Wardle-Konigl model ac- 
cretion disks. Our polarization maps depend on just 2 adjustable parameters, the 
grain axis ratio, a/b, and the mean ratio, £, of the energies in the random and ordered 
magnetic field components. 

2. We estimated the values of a/b and £ by performing a x 2 fit to the observed 100 /xm 
polarizations (H90, H93). Our models provide a reasonable fit to the observations if 
we assume that the la fluctuations in the observed polarizations due to turblence in 
the CND are > 0.6%. 

3. The observations are consistent with a continuum of models with different (a/b, £) val- 
ues. A model with a/b = 2/3 and £ = 1, the parameter values implied by independent 
constraints, falls close to the minimum of x 2 P er degree of freedom. For this model, 
X 2 /v = 0.97 and the mean error in the predicted polarizations is —0.2%. 

4. Our calculations imply that grains in the CND must be aligned with an efficiency 
such that (cos 2 (3) ~ 0.5, nearly independent of position. The degree and uniformity 
of alignment required by the observations are natural predictions of the ambipolar 
mechanism but require "fine tuning" of the parameters in other alignment theories. 

We thank Roger Hildebrand, Arieh Konigl, and the anonymous referee for helpful com- 
ments on earlier versions of the paper. This work was partially supported by NASA grant 
NAGW-3001. 
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Fig. 1. — Contours of x 2 / u — 1 (solid lines), 2 (short dash), and 3 (long dash) determined 
by fitting the grain axis ratio, a/6, and turbulence parameter, £. The filled circle marks the 
point (a/6, £) = (2/3,1) (see text). 
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—R.A. Offset (arc sec) 



Fig. 2. — Comparison between the far-infrared polarizations observed toward the CND (solid 
lines) and the predictions of a model with (a/b, £) = (2/3, 1) (dashed lines). Also shown are 
contours of the 100 /xm intensity predicted by our model, with contours at 20%, 40%, 60%, 
and 80% of the peak intensity. Offsets are measured from a = 17 h 42 m 29"4, 8 = -28°59 / 19" 
(1950). 



